Magnetic quantization of electronic states in d— wave superconductors 
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We derive a general quasiclassical approach for long-range magnetic-field quantization effects 
in superconductors. The method is applied to superclean d-wave superconductors in the mixed 
state. We study the delocalized states with energies e 3> Ao-y/ H/H C 2- We find that the energy 
spectrum consists of narrow energy bands whose centers are located at the Landau levels calculated 
in absence of the vortex potential. We show that transitions between the states belonging to the 
different Landau levels give rise to resonances in the a.c. quasiparticle conductivity and in the a.c. 
vortex friction. 



I. INTRODUCTION 



The unusual behavior of thermodynamic and trans- 
port properties of d-w&ve superconductors as functions of 
magnetic field is being a subject of extensive experimen- 
tal and theoretical studies. This behavior is attributed 
to nontriyial energy dependence of the electronic density 
of statesErta, and to specific kinetic processes which are 
very sensitive to fine details of electronic states brought 
about by the presence of vorticeaTQ. There exists, how- 
ever, a conceptual controversy about the structure of 
electronic states in <i-wave superconductors in the mixed 
state. One of the views is that the states below the maxi- 
mum gap Ao have a discrete spectrum due to Andreev l 
flections; some states are localized within vortex, cores 
while others are quantized at longer distanceal3ll3 as a 
particle moving along a curved trajectory in a magnetic 
field hits the gap for a current momentum direction such 
that e = A p . Other authors advocate that, instead of the 
magnetic quantization, energy bands should appear in. a 
periodic vortex potential due to the vortex latticeE3~E£l. 

In the present paper, we develop a general quasiclas- 
sical approach for calculating the long-range magnetic- 
field quantization effects in superconductors in the regime 
where the wave-length of electrons is much shorter than 
the coherence length pp£, S> 1. The proposed method 
is applied to superclean d-wave superconductors in the 
mixed state in the low field limit, H <C H c i- We demon- 
strate that quantization effects are in fact a compromise 
between the two abovementioned extremes. In the first 
part of the paper (Sections 0-0), we show that the 
influence of a magnetic field on delocalized excitations 
in a superconductor cannot be reduced to simply an ac- 
tion of an effective vortex lattice potential. The effect 
of magnetic field is rather two-fold: (i) It creates vor- 
tices and thus provides a periodic potential for electronic 
excitations, (ii) It also affects the long range motion of 
quasiparticles in a manner similar to that in the normal 
state. The latter long range effects are less pronounced 
for low energy excitations. The spectrum of excitations 



with energies e > /S.$yJ H / H C 2, however, is mostly deter- 
mined by the long range motion and exhibits magnetic 
quantization. 

W e study the delocalized states with energies e ^> 
Ao yj H/H C 2 and calculate their energy spectrum. We 
find that the spectrum indeed consists of energy bands 
as it should be in a periodic potential. However, in the 
quasiclassical limit, the bands are rather narrow; their 
centers are located at the Landau levels calculated in 
Refs. 

In the second part, Sections |v| - IVTl we consider effects 
of the energy spectrum on the vortex dynamics and on 
the quasiparticle conductivity. We show that both the 
vortex friction for oscillating vortices and the a.c. quasi- 
particle conductivity for fixed (pinned) vortices display 
resonances at transitions between the states belonging to 
different Landau levels. 



II. LONG-RANGE EFFECTS OF THE 
MAGNETIC FIELD 

We start with the conventional Bogoliubov-de Gennes 
equations 



2mA p w = 2meu . 



v — 2mA p u = —2mev (1) 



where p = — iV is the canonical momentum opera- 
tor. Equations (|l]) have the particle-hole symmetry such 
that u — > v*,v —* —u under complex conjugation and 
e — > — e. For a vortex array, the order-parameter phase 
is a multiple-valued function defined through 
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curlVx = 2tt5{v — n) . 



As a result 



1 



z x (r - r-j) 



(2) 



such that, on average, Vx ~ eHr/c for large r. 

Consider a quasiparticle in a magnetic field in the pres- 
ence of the vortex lattice for energies ranging from the 
above the gap to infinity. If the particle mean free path 
is longer than the Larmor radius, i.e., lo c t 3> 1 where uj c 
is the cyclotron frequency, such particle can travel away 
from each vortex up to distances of the order of the Lar- 
mor radius Tl = vp /lu c . This brings new features to Eqs. 
(|l|). Assume for a moment that A = 0. The wave func- 
tion u describes then a particle with the kinetic momen- 
tum P + = p (e/c)A and an energy e = P 2 i _/2m — Ep 
while v describes a hole with the kinetic momentum 
P = p+(e/c)A and an energy e = Ep — V'\j2ra. A par- 
ticle and a hole which start propagation from the same 
point will then move in different directions and along 
different trajectories which transform one into another 
under the transformation H — > — H. For a finite order 
parameter, the wave function is a linear combination of 
a particle and a hole. It is not convenient, however, to 
use such a combination at distances where the trajecto- 
ries of a particle and a hole go far apart, i.e., when the 
vector potential is no longer small compared to the Fermi 
momentum pp. 

Eq. (|l|) shows that the phase of u differs from that 
of v by the order parameter phase x- To construct a 
proper basis, one needs to bring the phases of u and v in 
correspondence with each other. We note that the usual 
transformation 



is not convenient when considering a particle which can 
move at distances much larger than the size of one unit 
cell. The problem is that the new functions u and v 
have extra phase factors ±x/2 as compared to the initial 
functions u and v, respectively. These phases increase 
with the distance resulting thus in a shift in the action 
A — ► A±x/2. The latter and is equivalent to a shift in the 
momentum p = VA — ► p± Vx/2. This transformation is 
not dangerous if the particle is bound to distances of the 
order of one intervortex distance because the phase gra- 
dient is limited |Vx| <C pf- However, for a vortex array, 
the phase gradient increases with distance and can reach 
values comparable with pp. It means that components 
of the new momentum can not be integrals of motion 
(i.e., they change along the trajectory) even in absence 
of the vortex potential associated with the superconduct- 
ing velocity and spatial variations of the order parameter 
magnitude. 

To avoid these complications we use another transfor- 
mation which also removes the coordinate dependence of 
the order parameter phase. The results, of course, should 
be independent of the choice of the transformation due 
to the gauge invariance. Following Refs. [ll] and [li] we 
put in Eq. (0) 



U = u, v = exp (— ix) v ■ (3) 
This is a single-valued transformation. We obtain 

u + 2me" !X Ap, v = 2meu , (4) 

v — 2me JX A^,, u = —2mev (5) 



P 2 -n 2 

1 Pf 



2mv, 



P% 



where P + = p — - A is the operator of the particle kinetic 
momentum, and 

P'+ = p Vx/2 = P+ mv s . 

The superconducting velocity is 

„ 2e , 

zmv s = Vx A . 

c 

In Eqs. (Q), (^|) we use that, for a general pairing sym- 
metry, Ap/ (x uv* depends actually on p' = (p u + p v )/2 
where p Ujt , are the canonical momentum operators which 
act on the Bogoliubov wave functions u and v, respec- 
tively. The term — Vx/2 appears in the order parameter 
together with the canonical momentum p because only 
one half of the momentum operator in A p / acts on each 
of the wave functions Morv. 

The transformation of Eq. (Q) is "w-like" and brings 
the phase of v in correspondence with the phase of u. 
The resulting equations are not symmetric with respect 
to u and v. the term v s is present in the second equa- 
tion together with P while it does not appear in the first 
equation. Let us perform one more transformation 



JXv/2 



(6) 



where Vx^ = 2mv s such that 

curlVx^ = V 2vr(5(r - rA - — H 

c 

i 

and Xv — X — Xa where 

XA = - fpxrVr'. (7) 

The "phase" Xv is not single valued within each unit cell, 
it depends on the particular path of integration. How- 
ever, it is single valued on average, i.e., on a scale much 
larger than the intervortex distance since 

curlVx^d 2 ^ = . 

It also implies that Xv does not have large terms increas- 
ing with distance. The transformation Eq. (JgJ) is thus 
not dangerous. The total transformation Eqs. (|^J^) has 
the form 



2 



u = exp (ix/2 - Ixa/2) U , 
v = exp (-ix/2 - ixa/2) V . 

With this transformation we finally obtain 



(8) 



P + — mv. 



|p + +mv s j -p% 
where 



U + 2mAp V = 2meU , 



V - 2mA p U = -2meV (9) 



e~ ix A 



p-(e/c)A — e2X ^-(e/c)A ■ 



As distinct from Eq. (Q) , a particle and a hole determined 
by Eq. (||) move along the same trajectory though, of 
course, in different directions. 

One can transform these equations further by putting 



# = 
where 



exp i I p • dr 



(10) 



(11) 



If divA = we have 



P+ (— iV + mv 5 ) 4>\ + mAp + 2 = rae0i , 
P+ (— tV — mv s ) 02 — "iAp + 0i = — me02 ■ 



(12) 



Another equation can be obtained using the transfor- 
mation 

u = e lx e- lx " ,2 U = exp (i X /2 + ixa/2) U , 

v = e - tx " /2 V = exp {-ix/2 + ixA/2) V . (13) 

We get 



P% 



P% 



U + 2mA p _V = 2meU, 

V - 2mA p _ U = -2meV (14) 



where P = p + (e/c) A is the "hole" kinetic momen- 
tum. The transformation Eq. ( |l3| ) is "w-like", it brings 
the phase of u in correspondence with that of v. Using 
Eq. (|l0|) we can transform Eq. (|l4|) to its quasiclassical 
version which is Eq. (|l2) where P + is substituted with 
P under the condition P_| 2 = p 2 F . Eq. (|l2] ) and its v- 
like analogue possess the particle-hole symmetry. Under 
transformation 



-e; 01 



they go one into another. Moreover, each set of equa- 
tions has the particle-hole symmetry separately for a 
given position on the trajectory if the kinetic momenta 



P± = p =F (e/c) A are reversed for a fixed position of the 
particle. Due to Eq. (11) p— {e/c) A — (gcosa, gsina), 
where a is the local direction of the momentum. The 
reversal corresponds to a — ► it + a. 

We take the z axis along the magnetic field and define 
the quasiclassical particle-like trajectory in Eq. (jl^) by 



dx 
dy 



Px - (e/c)A x 
Py ~ (e/c) A y 



(15) 



When the magnetic filed penetration length is much 
longer than the distance between vortices, S> ao> the 
magnetic field can be considered homogeneous. With A 
taken in the Landau gauge 

A = {-Hy, 0, 0) (16) 

the trajectory is a circle 

{x - x Q ) 2 + {y + c Px /eH) 2 = {p^c/eH) 2 (17) 

where p\ = p 2 F — p 2 . The local direction of the kinetic 
momentum is p x + eHy/c = p± sin a, p y = p± cos a. The 
distance along the trajectory is ds = r^da where the 
Larmor radius is Tl = p±/mu c . 

Eq. ([l2] ) has a simple physical meaning. It is the 
quasiclassical version of the usual Bogoliubov-de Gennes 
equation for vortex state modified to take into account 
long range effects of magnetic field. Eq. (02) c an be writ- 
ten in terms of the particle trajectory Eq. (|15[). We have 
fromEq. © 



'"j ( + mv t j 1 ' i -:■ A (n i < > 2 < • '1 . 



v± y-i-Q^ - mv tj 02 - A {a) 0i = -e^2 ■ (18) 

Here v± = px/m, and Vt is the projection of v s on the lo- 
cal direction of the trajectory. A(a) and v t are functions 
of coordinates x(s), y{s), and of the angle a{s) taken 
at the trajectory. Eqs. ([l8]) look exactly as the usual 
Bogoliubov-de Gennes equations. 



III. ELECTRONIC STATES IN ZERO LATTICE 
POTENTIAL 

For d-wetve superconductors, we take the order param- 
eter in the form A p = A {2p x p y ) / {p 2 + p 2 ) so that 

Ap-(e/c)A = Aosin(2a). Consider first the limit v s = 
and Aq = const. Eqs. (Es) become 



-iuj c — h A o sin(2a)0 2 = £01 • 

aa 

d(j) 

iu c — h A sin(2a)0i = e0 2 ■ 

aa 



With 



3 



4> = Cexp[i/(a)] 



we obtain 



/(«) = ± 



da 



e 2 - A 2 sin 2 (2a). 



The quantization rule also includes the integral over 
the momentum p defined by Eqs. ([l(], [ll]). We have 



p dr ± </> — 



e 2 - A 2 sin 2 (2a) = 2nn . 



(19) 



The quasiclassical approximation holds for n ^S> 1 . The 
± signs distinguish between particles and holes. As it was 
already mentioned, a particle (with the plus sign in Eq. 
(|l9|)) and a hole (with the minus sign) move along the 
same trajectory but in the opposite directions. The phase 
Xv which was introduced in Eqs. (^], ||) gives a contribu- 
tion to the action of the order of 2ir because it is limited 
from above by an increment of the order of circulation 
around one vortex unit cell; it can thus be neglected for 
large n. 



A. Sub-gap states 

In the range |e| < Ao, the turning points correspond 
to vanishing of the square root at a = ±a € where 
sin(2a e ) = |e|/A . We have 



4 



— daJe 2 - Ag sin 2 (2a) = 2?rn 



(20) 



where n > 0. The first integral in Eq. (19) disappears 



because the turning points of the momentum p are not 
reached: the particle can not go far along the trajectory 
Eq. ( |l5| ) and remains localized on a given trajectory at 
distances s ~ rt(e/Ao) smaller than the Larmor radius 
Th- Note also that the contribution from Xv vanishes 
identically because the particle after being Andreev re- 
flected transforms into a hole which returns to the start- 
ing point along the same trajectory. Using the substitu- 
tion sin a: = (Ao/e) sin (2a) we find 



daJe 2 - A 2 sin 2 (2a) 



A 



E 



e 

Ao" 
f_ 

A 2 



K 



A 



where K (k) and E (k) are the full elliptic integrals of the 
first and second kind, respectively. Ap plying the Bohr- 
Sommerfeld quantization rule Eq. (|19| ) we obtain 



2A n 



E 



Ao 



1 



A 2 



K 



A,, 



27m. (21) 



These states are degenerate with the same degree as in 
the normal state: for each n, there are <!>/2<i>o = N v /2 



states for particles and N v /2 states for holes, where $ is 
the total magnetic flux through the superconductor, and 
N v is the total number of vortices. 

Consider e <C Aq. Expanding in small k 



7T / k 2 \ , , 7T / k 2 



we find from Eq. (|2] 



e n = ±v/4A w c n. (22) 
Eq. (|2^) agrees with the result of Refs. [l(| |ll|. 



B. Extended states 

If |e| > A , we get for the Landau gauge Eq. ( |l6| ) 
p x = const and 

j> p dr = j> p y dy = 2 J \Jp\- (j>x + eHy/c) 2 dy 
= ircp\/eH . 

The turning points y\ t 2 correspond to the values of Lar- 
mor radius where Px~\~zHy\p/c — The correspond- 
ing trajectory is a closed circle where a varies by 2ir. The 
second integral in Eq. ( |l9| ) gives 

r^-A5si„S(2«) = i£ E (^V (23) 

The quantization rule ( |l9|) yields 

±h.E(^\=u c n+^--E F . (24) 
7T \ e„ J 2m 

For an s-wave superconductor we get, in particular, 

,2 



±Je 2 n -A 2 =u Jc n+£--E F . 



2m 



(25) 



IV. EFFECTS OF THE PERIODIC POTENTIAL 

A. Bloch functions 

At low magnetic fields H <C H C 2, one can consider 
that the particle trajectory always passes far from cores. 
The oscillating part of the order parameter comes mostly 
from the superconducting velocity. The corresponding 
Doppler energy rj = p±Vt is of the order of Aq^/ H/ H C 2- 
This periodic potential can split the energy spectrum into 
bands. Eqs. (0, |lj) or the quasiclassical version Eq. jl^ ) 
are invariant under the magnetic translations by periods 
of the regular vortex lattice. Consider the particle-like 
equations (|) or (12). The particle- like opepa,tor of mag- 
netic translations in a homogeneous field isli3 



4 



T (R/) = exp 



-iRj (p 



(26) 



where p = — iV is the canonical momentum and R/ is a 
vector of the vortex lattice. Its zero-field version corre- 
sponds to a shift 

f (Rj) / (r) = exp [-zRp] / (r) = / (r - R) . 

The operator T (R/) commutes with the Hamiltonian be- 
cause v s and A are periodic in the vortex lattice and the 
commutator 



c 



, (p- -J 

V c 



0. 



Since P+ does not change under the action of the oper- 
ator Eq. (|2q), magnetic translations for functions 4> in 
Eq. ( |l2] ) are equivalent to usual translations T (Ri) in 
space for a fixed kinetic momentum of the particle. 

It is more convenient to consider magnetic translations 
in the symmetric gauge A = H x r/2. In this case, 



T (R) / (r) = exp (^--R [H x r] f (r - R) . 

For this gauge, the wave functions Eq. ( |l0| ) can be more 
conveniently written in a slightly different form 



* (p x ; r) = exp 



ieHxy/2c + ip x x ■ 



y 

1 I Pydy' 
'm 



(27) 



The extra phase factor exp[ieHxy /2c] is associated with 
our choice of the vector potential and allows to re- 
duce the problem to the Landau gau ge. The parti- 
cle trajectory takes the form of Eq. ([17]) with p y — 



p\ — {p x + eHy/c) . The function <p satisfies Eq. (18) 
If ao and &o are the unit cell vectors along x and y 
respectively, the magnetic translation operators for func- 
tions of Eq. (p7|) are 



T x (la )V n (p x ; x,y) = e^^fe; x,y) , (28) 

eHlbo 



Ty(lb )^ n (p x ; x,y) = ^ n {Px 



x,y). (29) 



When deriving these expressions we have used the peri- 
odicity of v s and the fact that the trajectory depends on 
y only through y + cp x /eH. The turning point yi is thus 
shifted by Wq when p x is shifted by —eHlbo/c. 

The functions Eq. ( ^7| ) can be used to construct two 
independent basis functions 



*Z (k x ,k y ;x,y) -Y^,«.: 
$n ( k x,k y ;x,y) = 

xf y ((2l + l)b )y n (k x ; x,y) 



i 

E« 

i 



e ik « 2l( ">T y (2lb Q )$ n (k x ; x,y) , (30) 

ik v (2l+l)bo 

(31) 



with even and odd translations, respectively. Starting 
from Eq. (|3C]) we replace p x with k x . The functions $ ± 
belong to the same energy. The wave vector k y has an 
arbitrary value, we shall establish it later. The generic 
translation is 2bo which is the size of the magnetic unit 
cell along the y axis. The magnetic unit cell contains 
two vortices because the superconducting magnetic flux 
quantum correspond to one half of the 2ir phase circu- 
lation of a single-particle wave function. The functions 
Eqs. @ Hp) have the Bloch form 

f^lao^ik^ky) = (±l) l e-^ lao ^> ± (k x ,k y ), (32) 



Ty^mbo^i^.ky 



-ik v 2mbo ^± 



^ ± (k X ,ky 



(33) 



We omit the coordinates x, y in the arguments of $ ± 
for brevity. The functions $ transform into each other 
under odd translations 

f y ((2m + l)6o) $ ± (fc„ k y ) = e- ik ^ 2m+1 ^^(k X) k y ). 

(34) 

Since the magnetic translation T y (lbo) commutes with 
the Hamiltonian, the energy is degenerate with respect 



to k y . This degeneracy is spurious, however. To see this, 



consider the transformations Eqs 
the transformed function in Eq. (I 




For 



is periodic in k x 
with the period 2ir/ao. This period corresponds to the 
shift of the center of orbit yo = ck x /eH by one size of 
the magnetic unit cell 2&o- Obviously, the transformation 
Eq. (|3^) should also have the same symmetry. For one 
magnetic unit cell, a shift by 2b (i.e., for m = 1) along 
the y axis should combine with one period along the x 
axis. The period in k y is ir/bo; it should thus correspond 
to the shift of the coordinate xq by qq. This fixes 



k y = eHxo/c . 



(35) 



The energy depends on the position of the trajectory 
within the vortex unit cell through the Doppler energy rj. 
The energy e(k x , k y ) has a band structure due to period- 
icity of 77; it is periodic with the periods eHb^/c = tt/clq 
and eHcio/c = ir/bo in k x and ky, respectively, which 
correspond to shifts of the center of orbit by one vortex 
unit cell vector. 



B. Spectrum 



Consider energies e ^ Ao ^JH/H C 2- Applying the qua- 
siclassical approximation to Eq. (|l8|) we find 



(36) 



= Cexp [±iA(s)} 



where the action is 

A(s) = f ^/(e-r 1 ) 2 -Alsm 2 (2a)— . (37) 

J si ^-L 



5 



The quasiclassical approximation is justified because the 
wave vector dA/ds ~ efvp is much larger than the in- 
verse characteristic scale 1/ao of variation of the potential 
77 for e Aoy^H/ H C 2. The function ?; = p±Vt is taken 
at the trajectory which is a part of a circle specified by 
the coordinates of its center xq and y — ~ c Px/ e H; they 
determine the position of the trajectory within the vortex 
unit cell. 

For energies A \J Hj H C 2 <§; e < Ao, quasiparticle 
trajectory is extended over distances of the order of 
Tl (e/Ao)- The quantization rule defines the energy 

f 2 J(e- V f-A 2 sin 2 (2a)— = nn . (38) 
J Sl v ±- 

Here Si and S2 are the turning points. Expanding in 
small 77 -C e we find 



!/2 



/a / \le 2 - Agsin 2 (2a) — 
'2/1 ^ 

f V2 V(x,y)e dy 

m I — - — = 7m. 

•Mft Je 2 - A 2 sin 2 (2a) p f 



(39) 



Here r?(a;,y) = (/c x + eHy/c)v sx + p v v sy while j/ x and y 2 
are the turning points which correspond to vanishing of 
the square root: k x + eHy^/c = p^sin(2a e ). The en- 
ergy e n is a function of k x and Xq which determine the 
location of the particle trajectory with respect to vor- 
tices. The energy is thus periodic in k x with the period 
eHb /c and in xq with a period ao when the center is 
shifted by one period of the vortex lattice. 

The 77 term under the second integral in Eq. ( p9| ) os- 
cillates rapidly over the range of integration and mostly 
averages out. The remaining contribution determines the 
variations of energy with k x and xq and can be estimated 
as follows. For example, variation of action for e <C Ao 
due to a change in energy Se is 

SA ~ (Se/v F )(e/A )r L ~ (eSe)/(A cu c ) . 

Variation of action due to a shift of the center of or- 
bit by a distance of the order of the lattice period is 
5A ~ (oo/wf)?? ~ 1- The corresponding energy variation 
is thus Se ~ Ao^ c /e. Since Xq is coupled to k y through 
Eq. (|35|) the energy can be written as 



e n (k x , k y ) = JAAqoJc [n + 770 (k x ,k y )] 



(40) 



where 770 ~ 1 can depend on energy. The energy Eq. (|40|) 
has a band structure; the bandwidth is of the order of the 
distance between the Landau levels. It is small as com- 
pared to the energy itself. It is clear that the spectrum for 
energies e > Ao can also be obtained from Eqs. (pl|), ( gj ) 



and Q2q) through the substitution 71 > 71 -\- TjQ ( k x , ky ) 

One can check that, for given k x and k y , the quasipar- 
ticle states with different principle quantum numbers in 
deed concentrate near the levels determined by Eq. (|2 



if e 3> Aoy/H/H C 2- This is because the contribution 
to the action from the oscillating potential picked up on 
the distance of the order of the size of the unit cell ao is 
PfVsCio/vf ~ 1- The discrete structure of the levels Eq. 
( p2| ) would be preserved if the contribution to the ac- 
tion from the oscillating potential changes by an amount 
much less than unity for transitions between the neigh- 
boring levels. For an energy e, the distance between the 
neighboring levels is Se ~ Aqoj c / e. This corresponds to a 
change in the length of the trajectory by 

vp Se vp 

0S e -T • 

lu c A e 

The variation in the length is much smaller that the in- 
tervortex distance ao ~ £ \J H C 2 / H if e ^> Ao \J Hj H C 2 , 
and the action changes by a quantity much less than 1. It 
shows that the distance between the levels with different 
n is indeed determined by Eq. (p2|). 

The situation changes for smaller energies e < 

Ao \J H/H C 2'. The centers of bands will deviate strongly 
from positions determined by Eq. ( ^2[) due to a con- 
siderable contribution from theperiodic vortex potential 
to the turning points in Eq. (|38|). Moreover, the appli- 
cability of the quasiclassical approximation, i.e., of Eq. 
( |38| ) itself is violated; the potential r\ is strong enough to 
cause large deformations of the energy spectrum. As was 
shown in Ref. ^ some states can even become effectively 
localized near the vortex cores. 



V. INDUCED TRANSITIONS BETWEEN THE 
LANDAU LEVELS 

Vortex motion induces transitions between the quasi- 
particle states. The transitio ns between low-energy core 
states with e -C Ao \J H/ H c % were considered in Ref. [| It 
was shown that the vortex core states determine the vor- 
tex response to d.c. and a.c. electric fields. For tempera- 
tures T C \J H I H C 2 <§; T extended states dominate. It was 
found in Ref. [| that the vortex response is determined by 
what was called "collective modes" which are associated 
with the electron states outside the vortex cores. In this 
Section we demonstrate that these collective modes are 
nothing but transitions between the electronic states Eq. 
( ppj ) specified by the same quasimomentum but by dif- 
ferent principal quantum numbers n. We start with not- 
ing-that the transition matrix elements arc proportional 



($„ (k,) V-Hi$ m (kj )) where the Hamiltonian Hi is 
composed of Ap and 77, while k is the quasimomentum. 
Hi is periodic with the period of the vortex lattice thus 
the transitions are possible between the quasimomenta 
which differ by vectors of the reciprocal lattice. Since 
the band energy is periodic in the quasimomenta with 
the periods of the reciprocal lattice, the energy difference 
for these transitions corresponds to the energy difference 
for states with the same quasimomentum but with dif- 
ferent quantum numbers n. For 770 <C n the transition 
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energy is just the distance between the Landau levels: 
5e n (k x ,k y ) = 6e n determined by Eqs. ( [2l"| , ^J) or (25). 
For low energies in a d-wave superconductor, one has 
de (k x , k y ) — 2AodJ c /e„ in accordance with Eq. ([22]). 
Consider transitions between the levels which are ex- 



cited by oscillating the vortices in more detail. We use the 
microscopic kinetic-equation approach which has been 
applied earlier for s-wave superconductors in Ref. [n]. 
The kinetic equations, for the distribution functions fi 
and f 2 have the formE2l 



e (v F • E) 5 _ 



1 f ^ +f t^E 

2 \ J ~ at J ~ dt 



- [v F x H] g_ } 



_VA* + /IVAj 




df (0) afi 

^ r + (v,.V)(. 9 _f 2)+5 ^ 



1 f f 9A; 

2 y- dp 



dp 



Vfi = J 



(41) 



and 



Here 



ff _ (v F • V) fx = 0. 



(42) 



MA) 



g R(A) fR(A) 
_f\R{A) 



-<] 



R(A) 



are the retarded (advanced) quasiclassical Green func- 
tions, and g- — [g R — g A ) /2. 

For an extended state with an energy e > A p , the par- 
ticle trajectory crosses many vortex unit cells at various 
distances from vortices. Since the distribution function 



fi is constant along the trajectory according to Eq. (42), 
it should be also independent of the impact parameter 
(i.e., of the distance from the trajectory to the vortex). 
We thus look for a distribution function f\ which is inde- 
pendent of coordinates. One can then omit the last term 
in the l.h.s. of Eq. (f4l|). Let us average Eq. (^TJ) over an 
area which contains many vortex unit cells but has a size 
small compared with the Larmor radius, ao <C r <C tl. 
Since r « ri the momentum p is still an integral of 
motion. We have (compare with Ref. |20|) 



So 



1. 



dp 



Jd A r 



So 



-Tr / dfrg- (X7H 
2 Jso 

" 1 d 2 rg-(v L -VH" 



-Tr 



So 



dc 



Here Tr is the trace in the Nambu space, So = $o/B is 
the area of the vortex unit cell, 



J 



1 



[(fr <<?->- &<?-»<?- 



h (/i) - (fi/i)) /- + (h (/-) - £ 



Using the identity 

±Tr f d 2 r [(VH) S-] 

J On 



7T [z X V_|_] 



-7T [z X VjJ 



dh dfi 
dp dt 



.9- d r 



So 



Jd z r 



So 



= 7T (V L • [Z X VjJ) 



gfW 
de 



(43) 



We shall concentrate on energies e 3> A^/H/H^. In 
the leading approximation 



9- 
f- 



V^ 2 ~ A 2 (a) 

A(a) 
^e 2 - A 2 (a) 



9[e 2 
8[e 2 



A 2 (a)] , 
A 2 (a)] . 



We have 



(fx) = (f l5 _> = 0; (f 1 /_) = (f 1 /l) = 0. 

For a d-wave superconductor also (/_) = ^/l^ = 0. 

In the collision integral, the main contribution for 
e ^> H / H C 2 comes from the delocalized states. In- 
deed, including contributions from the bound states in 
the cordis with energies E n (b) we would have 



Jdfr 



S 



So 



{g-)g- 



fi 



where b is the impact parameter. The first term in square 
brackets comes from the core states. Since r n ~ t and 
b ~ £,y H C 2/H, the core contribution is of the order 
of t -1 y H/H C 2- The delocalized states, however, give 
(e/Ao) t -1 which is much larger than the first term. Ne- 
glecting the core contribution we find 



J = 



(g-)g-fi- 



Let us put 

9f (o) 



fi 



de 



[([u x pj • z) 7o + (u • p x ) 7H ] (44) 



derived in Ref. 19 we find 



The functions 70, h satisfy the following set of equations 
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- 7h - V {a) 70 + 1 = 

-5 1- 7o - V" (a) 7h = 

oa 

which is derived from Eq. ([43|). Here 

(-iu> + (g_) /r)g- 



V(a) 



(45) 



(46) 



The general solution of Eqs. (|4J) can be obtainedH by 
putting W± = 7h ± «7o ■ We have 



dW± 
da 

The solution is 



=F zW± - V (a) W± ± « = . 



W± = 
where 



Jo 



(47) 



F(a)= V (a') da' . 



The constant C± is found from the condition of period- 
icity W (a) = W(a + Tr/2) 



exp [F (n/2)] J^ 2 exp [+ia - F (a)] da 



C± 



1 - exp [±Z7r/2 + F (V/2)] 
In the limit t — + 00, the solution 

7h = {W+ + WJ) /2; 70 = (W+ - WJ) /2i 
has poles when 



F(7r/2) = ^-(l + 2M) 



(48) 



(49) 



where M is an integer. If e > Ao, one obtains resonances 
at 



da = 27r(l + 2M) 



The lowest frequency M = exactly corresponds to the 
condition 

u> = (de n /dn) 

where de n /dn is the distance between the Landau levels 
determined by Eq. fll9j). 

Note that, for an s-wave superconductor, Eqs. (|45|) 
has the form 



7h + V"7o = 1 
7o - Vjn = 



(50) 



where 



V(a) 



6 [ £ 2 - Aq] 



1 



J= 6 [e 2 - A*] 



One has from Eq. (|50j) 

7h 



1 



l + V 2 



70 



l + V 2 



The resonances appear when lo c t 3> 1; the poles corre- 
spond to V = ±i so that 



u; = W, 



^e 2 - Aq _ de r , 



where e„ is determined by Eq. (p5|). 

A. Low energies 

For energies e < A , the resonance condition Eq. ( |49| ) 
is not just the distance between the Landau levels deter- 
mined by Eq. (f2l|). One has from Eq. ( |2(i| ) 



=da = 



(in 7- Qe V? - A 2 (a) 

where A (a e ) = e. At the same time, Eq. (|i^) gives the 
lowest resonant frequency 



— N / - 



da = 2-7T 



where IV is the number of gap nodes ( N = 4 for a d- 
wave superconductor) . We see that the resonance occurs 
at 



Nlu = 2 



dc ra 

dn 



(51) 



When the vortex oscillates, all N nodes participate in 
exciting quasiparticles which accounts for the factor N 
111 the l.h.s. of Eq. ©. This is similar to the process of 
multi-photon absorption. The factor 2 in the r.h.s. is ex- 
plained by noting that states with momentum directions 
a and a + n are simultaneously excited. 

Solution of Eqs. ©H) for A ^H/H c2 < e < A 
was obtained in Ref.TJ. For the main re gion of angles, 
\a\ > a e — e/2A . According to Eqs. @^|), it is 



70 = A cos a + B sin a 
7h = 1 — A sin a + B cos a 



with 



A = 



e sinh A 
2sinh 2 A+ L 
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' sinh A 



2 sinh A + 1 



(52) 



(53) 
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Here we use that F (tt/2 



2A — F (a) where 



VI. QUASIPARTICLE CONDUCTIVITY 



A = F{a € ) ;F(?r/2) = 2A. 



One has 



A = 



l/r e// 7T [e| 



4A 



where \/r e ff = |e| /Aot since (#_) = ]e| /Aq. Note that 
a r-approximation was used in Ref. [| for the collision 
integral. To get the present expression for A from that 
obtained in Ref. ^ one has to replace 1/r with l/r e ff. 
For r — > oo, the response Eqs. ([32|), ( pq ) has poles at 
iX = (2M + 1)tt/4, i.e., for 



w = (2M + l)S (e); B (e) = A ^ c / |e| 



(54) 



We have for M = 



_ 1 de n 



where e„ is determined by Eq. (p2]). This condition 
agrees withEq. @. 

These resonances were first predicted in Ref. |9|. Note 
the different numerical factor in Eq. (^4j) as compared to 
Ref. [)| this is because a simplified version of V (a) has 
been used in Ref. ^. The main effect of resonances is that 
vortices experience a considerable friction force Eq. (pi 
even in a superclean case wr e // 3> 1. 



Consider the a.c. quasiparticle conductivity which can 
be observed if vortices are pinned. The distribution func- 
tion can be found from Eqs. (41, fl2|). We are looking 
again for the distribution function f\ which is indepen- 
dent of coordinates. One has 

e (v F ■ E) 9--g^ + (vf • V) (,g_/ 2 ) + ,g_ — 



; [v F x H] g_ - 1 (/_ VA£ + /lVA p 



dp 



J. 



We omit the time derivatives of A because vortices do 
not move. After averaging over the vortex lattice we get 



r i 5fl 9f i f ^ 
n[zXV±] -d^--dtJ So 9 - dr 

df (o) 



Jd 2 r 



So 



df (0) f 



For the distribution function in the form 



h = 



gfW 

<9e 



[(E-p ± )7o-([ExpJ-z)7H] (58) 



B. Vortex friction 



we obtain 



A vortex moving with a velocity V£ experiences a force 
from the environment 



-Dv L - D'[v L x z] 



(55) 



According to Ref. |l9|, the delocalized states contribute to 
the friction constant 



d/(°) de 



Ddel = 7VN{ L 7 ° de 2 



(56) 



where (• • •) is an average over da. The factor D 1 is de- 
termined by the same expression where 70 is replaced 
with 7h- 

The presence of resonances makes the dissipative con- 
stant Z?dci finite even in the superclean limit lo c t — * 00. 
Indeed, for an s-wave case, 



70 



^- [<J(w -E) + S(uj + E)] 



where E 
comes 



0J c y/l — Ap/e 2 . The friction constant be- 

W 2 K 2 # (0) M 



D de i = TT NAq 



(1 - ujy^f 2 



de 



(57) 



where eo = / v/l — ^V^c- more detailed discussion 
of the resonant vortex friction for a d-wave superconduc- 
tor at low temperatures can be found in Ref. |. 



7H - K (a) 70 = g- 

oa muir 



9th 



9a 

The solution is 



+ 7o - V (a) 7h = . 



W± 



ic: 



moj c J 



9-e 



Tta >-F( a ') da , + c 



(59) 



±ia+F(a) 



(60) 



The periodicity condition W (0) = W(ir/2) gives for 
e< A 



C*4 



1 



to [— iu> + (g_) /7 



sinh A / cosh A 
±i 1 



cosh 2 A 



cosh2A 



(61) 



where A = F'(a e ). 

The quasiparticle current is 



j(«0=-„(0)e Jvpg^de^- 
= ^ p, E + 4 ?!,) [Exz] 
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where 



(gp) _ ^(0)e 



y O,H 



v±p±g-7o,H- 



de. 



2 J '"'"4tt <9e 

Calculating the integral over da in Eq. (|62|) we find 



<7oS- 
(7h5- 



4mu>r X 2 



X - 



tanh A 
tanh 2 A + 1 



7re (g-) 2 tanh 2 A 
4muj c A 2 tanh 2 A + 1 ' 



(62) 

(63) 
(64) 



and 



A = ^ (5-) + (3-) /t] ^ c 1 • 

Since A is independent of the momentum directions, the 
quasiparticle conductivity becomes 



Jip) 



Ne / (5-7o) 



0f W de 
de ~2 



(65) 



and the same expression for Ojj where (g-7o) is re- 
placed with (<7_7h) • 

Consider first the superclean limit u; c T e // ^> T/T c such 
that ReA < 1. For wr e // > 1, where r e// - (T c /T) r the 
response Eq. ( |63| ) has resonances at iA = (2M + 1)tt/4 
which is again the condition of Eq. (54): 



(5-7o) = 



4e 



,t/,), ; - c (2Af + 1 



The dissipative part of the quasiparticle conductivity be- 
comes 



It is 



w) 2JVe 2 w 2 Ao 



TTtnT 



cosh 



M=0 



A w c (2A/ + 1) 



2TU 



Recr 



(<?p) 



2iVeV 



for u ^> Eg where E g = Aoo-> c /2T, but decreases expo- 
nentially for smaller u ^ E g . The dissipative part for 
uj <C Eg is mostly due to r. Since A -C 1 one has in this 
limit 



Rcct, 



(?p) _ 
o — 



irNe 2 
3mu> c 



df<®de. , 2 . 



5tt 7 T 4 



24mw 2 rA^ 



On the moderately clean side, such that w c r e // <C 
T/T c one has ReA ^> 1. The conductivity has a Drude 
form 



(qp) _ 
a Q — 



Ne 2 f 0f(°> de 



{9-) 



de 2 [-iuj+(g_)/T] 



The Hall conductivity does not contain contributions 
from poles because the resonances in Eq. (|64|) with 
M > cancel those with M < 0. For lo <C (<7_) /r 
one has 



(7Hff-) Q 



tanh 2 - 



m (tanh 2 w + l)w 
where w (e) = it (g~) 2 /4u; c t. The conductivity is 
Ne 2 r f°° Of (°) tanh 2 to de 



(gp) _ 

CT H — 



<9e tanh 2 u> + 1 w 



If T <C Aoy/uj c T one has e ~ T and w -C 1. In this 
limit 



Mp) 



nNe 2 
Amu, 



9f(°) , x2 , ttWT 2 
(9-) de = 



c JO 



12mcij c A 2 ' 



If T 3> Ao^/w c t and w c r <C 1 the integral is determined 
by w ~ 1 and e ~ Ao-^/ZZvr <C T. We have 



0.39Ne 2 TA ^/UJ~T 
Tm 



VII. CONCLUSIONS 

We discuss and analyze the "Landau levels" vs "energy 
bands" opposition concerning the structure of the excita- 
tion spectrum in the mixed state of superconductors, and 
in particular, of d-wave superconductors. We find that 
the actual picture of quantization is an interplay between 
the two limiting images of the energy spectrum. Our 
analysis shows that the influence of the magnetic field on 
delocalized excitations in a superconductor can not be 
reduced to a mere action of the effective vortex lattice 
potential. In fact, magnetic field has a two- fold effect: 
On one hand, it creates vortices and thus provides a pe- 
riodic potential for excitations, on the other hand, it also 
affects a long range motion of quasiparticles in a manner 
similar to that in normal metals. For low energy excita- 
tions, the long range effects are less pro nounced . How- 
ever, excitations with energies e > Aoy/ H/ H C 2 mostly 
show the long range quantization. The energy spectrum 
consists of "Landau levels" which are split into bands by 
the periodic vortex potential. In the quasiclassical ap- 
proximation pp£ 3> 1, the bandwidth is of the order of 
the distance between the Landau levels; it is small com- 
pared to the energy itself. 

An a.c. electric field induces transitions between the 
states belonging to different Landau levels. Using the 
microscopic kinetic equations we demonstrate that these 
transitions can be seen as an increase in the vortex fric- 
tion and/or in the quasiparticle conductivity due to a 
resonant absorption at frequencies corresponding to the 
energy differences between the Landau levels. 
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